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SUMMARY:  Damage  evolution  in  a  particulate  composite  around  a  notch  tip  was  studied. 
The  composite  had  hard  and  stiff  particles  embedded  in  a  soft  and  weak  matrix  like  a  rubber 
material.  The  major  damage  modes  were  the  interface  debonding  called  dewetting  and  the 
resulting  matrix  cracking  because  the  particles  were  much  stronger  than  the  matrix.  A 
numerical  modeling  and  simulation  of  such  damage  was  conducted  using  the  micro/macro¬ 
approach.  This  technique  combined  micro-level  analysis  and  macro-level  analysis.  Damage 
was  described  at  the  micro-level  using  a  damage  theory.  Damage  initiation  and  growth  at  a 
circular  notch  tip  were  predicted  from  the  numerical  study,  and  their  results  were  compared -le^ 
experimental  data.  Both  results  compared  very  well. 


KEYWORDS:  damage,  particulate  composite,  micromechanics,  finite  element  analysis, 
micro/macro-approach 


INTRODUCTION 

Particle-reinforced  composite  materials  have  been  ‘used’ .  widgly.  in  various  engineering 
applications  ranging  from  commercial  applications  to  military  applications.  Their  damage 
modes  and  mechanisms  depend  on  properties  of  the  reinforced  particles,  the  binder  matrix, 
and  their  interface.  In  general,  the  damage  modes  in  a  particulate  composite  are  matrix 
damage  (micro-cracking),  interface  debonding,  and  particle  cracking.  Relative  strength  and 
stiffness  of  the  constituent  materials  and  their  interface  including  initial  defects  will  govern 
the  damage  modes.  The  composite  material  under  study  has  hard  and  strong  particles 
embedded  in  a  very  soft  and  weak  rubber-like  matrix  material.  As  a  result,  the  major  damage 
modes  are  interface  debonding  and  matrix  cracking  [1].  In  particular,  the  interface  debonding 
called  dewetting  is  the  primary  damage  mode. 

Previous  research  on  these  kind  of  composite  materials  has  focused  on  material  degradation 
based  on  damage  evolution.  Those  studies  attempted  to  develop  macro-level  constitutive 
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models  incorporating  progressive  damages  [2,3].  Other  studies  considered  the  progressive 
damage  from  notch  tips  [4-7].  Most  of  the  previous  studies  did  not  consider  dewetting 
explicitly.  It  was  treated  as  a  smeared  damage  in  the  composite  [2-4]  or  in  the  matrix 

[5-7]- 

The  objective  of  the  paper  is  to  model  the  progressive  interface  damage  starting  from  a  notch 
tip  using  a  micro/macro-mechanical  approach.  The  interface  damage  (dewetting)  is  considered 
explicitly  at  the  micro-level  analysis  of  the  micro/macro-approach.  The  subsequent  section 
describes  the  micro/macro-approach  including  the  macro-level  analysis,  micro-level  analysis, 
and  a  damage  theory.  Numerical  results  and  discussion  are  presented  next  followed  by 
conclusions. 


MICRO/MACRO-APPROACH 

The  modeling  and  simulation  technique  used  in  this  study  is  called  the  micro/macro-approach 
[5-10].  The  approach  is  similar  to  the  local-global  approach  in  structural  analysis.  The 
micro/macro-approach  utilizes  the  two  levels  of  analyses:  micro-analysis  and  macro-analysis. 
The  micro-analysis  is  performed  at  the  constituent  material  level  while  the  macro-analysis  is 
undertaken  at  the  composite  structure  level.  The  two  analyses  are  conducted  in  tandem  as 
damage  initiates  and  evolves  in  composite  structures.  Damage  is  described  at  the  micro-level 
in  terms  of  the  constituent  materials.  Damage  modes  for  a  particulate  composite  are  described 
in  terms  of  particle  cracking,  matrix  micro^cracking,  and  interface  debonding.  The  interaction 
between  the  two  levels  of  analyses  is  shown  in  Fig.  1  and  the  description  of  each  analysis  is 
given  below. 


MACRO-ANALYSIS 
=>  FEA  of  Composite 
Structures 


Composite  Strains 
- - ► 

< - - - 

Effective  Property 


MICRO-ANALYSIS 
=>  3-D  Unit-Cell  Model 


Micro-Strains 


Damage  State 
Degraded  Property 


r 


DAMAGE  THEORY 


Fig.  1:  Interaction  between  micro-analysis  and  macro-analysis 

Macro-Analysis 

The  macro-analysis  is  conducted  at  the  composite  structure  (specimen)  level.  Thus,  the 
geometry  and  dimensions  of  the  structure  are  at  the  macro-scale.  For  this  level  of  analysis,  any 
structural  analysis  technique  can  be  adopted.  However,  in  order  to  solve  general  cases  of 
composite  structures  under  diverse  loading  and  boundary  conditions,  numerical  solution^ 
techniques  are  suitable  such  as  the  finite  element  method  [11-13],  boundary  element  method 
[14,15],  finite  difference  method  [16],  and  meshless  method  [17].  Some  combinations  of  these 
techniques  may  be  used  to  take  advantage  of  respective  technique*?  Because  the  finite  element 
method  has  been  most  popular  for  structural  analysis,  it  is  used  for  the  present  study. 


The  composite  structure  under  study  is  discretized  into  a  finite  element  mesh  with  proper 
structural  elements  such  as  plate/shell,  beam,  2-D  continuum,  and/or  3-D  continuum  elements. 
In  general,  isoparametric  types  of  elements  are  used  with  the  Gauss  quadrature  rule.  These 
elements  require  material  properties  at  the  Gauss  quadrature  points.  As  a  result,  the  interaction 
between  the  micro-analysis  and  the  macro-level  finite  element  analysis  is  performed  at  the 
Gauss  quadrature  points.  The  micro-analysis  provides  the  Gauss  quadrature  points  of  the 
macro-analysis  with  effective  composite  material^ properties  including  the  damage  state.  The 
macro-analysis  is  undertaken  to  compute  deformation  and  stress/strain  of  the  composite 
structure.  These  data  are  transferred  to  the  micro-level  analysis. 

The  present  composite  structure  (specimen)  to  be  studied  is  a  square  plate  (76.2mm  x  76.2 
mm)  with  a  through-the-thickness  hole  with  diameter  of  6.35mm  at  the  center.  Because  of 
symmetry,  a  quarter  of  the  plate  is  modeled.  The  specimen  is  subjected  to  a  uniaxial  tensile 
load.  A  uniform  displacement  is  applied  to  the  edge  of  the  specimen  so  that  the  numerical 
results  can  be  compared/ttf  the  experimental  data. 

Micro-Analysis 

Because  the  micro-analysis  is  conducted  repeatedly  at  each  Gauss  quadrature  point  of 
,  respective  finite  element  with  iterative  processes  damage  progresses,  the  analysis  should  be 
computationally  efficient.  To  this  end,  a  simplified  analytical  model  is  used  for  the  micro¬ 
level  analysis.  One  of  them  is  a  unit-cell  model  [18-22], 

The  micro-analysis  has  two  major  functions.  One  is  to  compute  effective  material  properties 
of  the  composite  based  on  the  progressive  damage  (af)the  constituent  materials.  These  data  are 
used  for  the  subsequent  macro-analysis.  The  other  function  is  to  decompose  the  macro-level 
stress/strain  obtained  from  the  previous  macro-analysis  into  micro-level  stress/strain.  The 
micro-stress/strain  is  used  to  determine  the  damage  state  at  the  constituent  material  level.  The 
continuum  damage  theory  is  utilized  to  describe  the  damage  state  as  well  as  the  degraded 
material  properties  of  each  constituent  material.  The  micro-analysis  utilizes  the  degraded 
material  properties  to  calculate  the  updated  material  properties  of  the  damaged  composite. 

The  micro-model  is  shown  in  Fig.  2.  It  is  a  unit-cell  consisting  of  eight  sub-cells.  The  unit  cell 
has  dimensions  of  unity  in  every  direction.  One  sub-cell  (sub-cell  fl  in  Fig.  2)  is  the  particle 
sub-cell  and  it  is  a  cubepFside  length  of  (Fp)1/3  where  VP  is  the  particle  volume  fraction  of  the 
composite.  Then,  sizes  of  the  rest  of^ub  -cells  are  determined  accordingly.  In  order  to  model 
the  particle/matrix  interface  debonding,  linear  springs  are  introduced  at  the  interfaces  of  the 
particle  sub-cell  and  three  neighboring  sub-cells.  (The  double  arrows  in  Fig.  2  indicate  linear 
springs  at  the  interfaces)  Before  interface  failure,  the  linear  spring  is  very  stiff.  However,  as 
interface  damage  called  dewetting  evolves  in  the  model,  the  spring  constant  decreases 
accordingly.  The  amount  of  reduction  in  the  spring  stiffness  is  dependent  upon  the  damage 
state.  This  will  be  discussed  in  the  next  section. 


For  mathematical  simplicity  that  is  an  important  part  of  the  success  of  the  micro/macro¬ 
approach  in  a  practical  sense,  stress  and  strain  are  assumed  to  be  uniform  within  each  sub-cell. 
Then,  stress  equilibrium  at  the  interface  of  any  two  neighboring  sub-cells  is  enforced. 


_  3  4  _  5  __6  _7  8 
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2  _  4  5  _  7  6  _  8 

C  22  —  ^22  ^  22  ^  22  ^*22  —  ^22 


(1) 

(2) 


(3) 


in  which  superscript  denotes  the  sub-cell  number  as  shown  in  Fig.  2  and  subscript  indicates 
the  stress  component.  These  equations  are  for  normal  components  of  stresses.  Similar 
equations  can  be  written  for  shear  components  but  they  are  omitted  here. 

The  deformations  of  sub-cells  are  assumed  to  satisfy  the  following  relationship: 


^>£n  +  C£  11  +  (^an  /  n  +  4>£n 

22  4i£22  +  ([p®22  '  ^2  ^  ^>£ 22  4>£ 22 

^£ 33  +  C£ 33  ({p® 33  I  lp£ 33  +  ^m£ 33 


=  ^>£  11  +  4i£11  =  ^Sll  +  4iSll 

=  ^/>£22  +  ^nS22  =  ^7£22  +  4i£22 
=  ^p£33  +  ^m£33  =  ^/>S33  +  4>£33 


(4) 

(5) 

(6) 


where 


/, = vr  (7) 

(8) 

Here,  Vp  is  the  particle  volume  fraction  of  the  composite.  In  Eqs^(4)-(6),  kx ,  k2 ,  and  kz  are 
spring  constants  in  the  respective  interface  as  shown  in  Fig.  2.  For  isotropic  damage,  these 
constants  are  assumed  to  be  the  same. 
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Other  necessary  mathematical  expressions  are  the  constitutive  equations  for  the  particle  and 
the  matrix  as  well  as  for  the  composite.  'Furthermore,  the  composite  stress  and  strain  are 
assumed  to  be 
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(9) 


(10) 


where  superimposed  bar  denotes  the  composite  (macro)  stress  and  strain,  and  V"  is  the 
volume  of  the  n*  sub-cell. 

Algebraic  manipulation  of  the  previous  equations  finally  yields  the  two  key  expressions  as 
given  below: 


(li) 

{<s}=[R2]{s'ff}  (12) 

in  which  is  the  effective  composite  material  property  matrix,  [V\  is  the  matrix  composed 
of  sub-cell  volume  fractions,  [E\  is  the  matrix  consisting  of  constituent  material  properties, 
and  [R2]  is  the  matrix  relating  the  sub-cell  strain  vector  consisting  of  particle  and  matrix 
strains  {s}to  the  unit-cell  strain  vector  (effective  composite  strain)  {seff}. 

Equation  (1 1)  is  used  to  compute  the  effective  material  properties  of  the  composite  material 
based  on  constituent  material  properties  and  their  geometric  information.  Any  degraded 
material  property  of  the  constituent  materials  associated  with  damage  evolution  is  used  to 
compute  the  effective  property.  Thus,  this  equation  provides  the  connection  from  the  micro¬ 
analysis  to  the  macro-analysis.  On  the  other  hand,  Eq.  (12)  provides  the  connection  from  the 
macro-analysis  to  the  micro-analysis.  That  is,  the  equation  decomposes  the  effective 
composite  strains  into  micro-level  strains  at  the  particle  and  matrix.  The  micro-strains  are 
used  for  the  damage  theory  as  described  below  •  to  determine  the  progressive  damage, 

and  they  are  also  used  to  compute  micro-stresses. 

Continuum  Damage  Theory 

Damage  at  the  micro-level  is  described  using  a  continuum  damage  theory.  In  general,  an 
anisotropic  damage  theory  can  be  applied  to  the  constituent  materials  because  the  damage 
state  may  be  anisotropic  depending  on  the  given  conditions.  An  anisotropic  damage  theory, 
however,  requires  many  material  parameters  most  of  which  are  not  easily  available.  On  the 
other  hand,  a  simpler  isotropic  damage  theory  may  be  used  at  the  micro-level  assuming  that 
each  constituent  material  behavior  is  isotropic  at  the  level  and  its  damage  state  is  also 
isotropic.  These  assumptions  are  realistic.  Therefore,  an  isotropic  damage  theory  is  used  in 
this  study. 

Because  the  interface  debonding  called  dewetting  is  the  major  damage  mode  in  the  present 
composite,  the  present  damage  theory  addresses  this  kind  of  damage.  It  is  assumed  that  the 
damage  growth  at  the  interface  is  related  to  the  change  of  the  volume  like 

d<j)=(  1  — 4>)  (13) 

where  <j>  is  the  damage  parameter  and  stl  is  file  strain  tensor  of  a  material.  This  equation  applies 
when  an  induced  strain  measure  exceeds  the  threshold  value  which  is  material  dependent.  The 
damage  expression  is  based  on  the  experimental  observation  which  showed  that  the  volume 
dilatation  increases  as  damage  grows. 

Once  the  damage  parameter  <|>  is  computed,  the  interface  spring  constant  is  adjusted  based  on 
the  amount  of  the  damage  using 


(14) 


in  which  k  is  the  reduced  interface  spring  constant  and  k0  is  the  original  spring  constant  before 
interface  damage.  The  function /is  assumed  as 

/♦)  =  10"“  (15) 

where  c  is  a  material  constant.  Another  kind  of  damage  theory  was  used  in  Re£  4  and  6. 

RESULTS  AND  DISCUSSION 

First  of  all,  the  stress-strain  behavior  of  a  rniiform  specimen  of  the  composite  is  plotted  in  Fig. 
3.  The  experimental  data  are  compared  ie  the  numerical  data  obtained  from  the  micro/macro¬ 
approach.  The  two  curves  compare  well  until  final  fracture  of  the  specimen.  Damage 
initiation  and  growth  is  shown  in  Fig.  4.  As  seen  in  the  figure,  damage  initiates  at  the  applied 
strain  level  of  a  little  less  than  0.05  m/m,  where  Fig.  3  shows  the  deviation  from  the  linear 
stress-strain  curve.  That  is,  damage  initiation  and  growth  result  in  nonlinear  stress-strain 
behavior  until  fracture.  The  maximum  value  of  damage  parameter  in  Fig.  4  is  considered  as 
the  saturation  of  damage. 


Fig.  3:  Stress-strain  curve  of  a  uniform  composite  specimen 


In  general,  damage  growth  in  a  uniform  composite  specimen  is  not  uniform  across  the 
specimen  because  of  non-uniform  micro  structures  including  nonfuniform  particle  distribution 
in  the  composite.  Thus,  aamage  varies  from  point  to  point  in  the  specimen.  However,  the 
analytical  model  cannot  implement  such  a  non-uniform  variation.  Instead,  only  the  average 
damage  in  the  composite  can  be  simulated.  As  a  result,  a  local  maximum  damage  value  in  the 
composite  specimen  may  be  greater  than  the  maximum  value  shown  in  Fig.  4.  In  the 
subsequent  study,  this  maximum  damage  value  is  used  as  the  damage  saturation  value  and  the 
percentage  of  damage  is  computed  based  on  that  maximum  damage  value. 

The  next  example  is  a  composite  specimen  with  a  through-the-thickness  circular  hole  at  the 
center  of  the  specimen.  The  specimen  geometry  was  provided  in  the  previous  section.  The 
specimen  is  subjected  to  a  uniform  displacement  along  two  opposite  edges. 

Figure  5  shows  the  plot  of  the  applied  strain  versus  the  induced  strain  at  the  notch  tip.  The 
applied  strain  is  computed  from  the  applied  displacement  divided  by  the  specimen  length. 
The  notch  tip  strain  is  the  component  along  the  loading  direction.  As  expected,  the  induced 
strain  at  the  notch  tip  is  three  times  greater  than  the  applied  strain  due  to  strain  concentration 
at  a  low  applied  strain  level.  As  damage  initiates  and  grows  at  the  notch  tip,  the  material  at 
the  tip  degrades.  As  a  result,  the  induced  strain  increases  greater  than  the  applied  strain.  That 
is,  the  slope  in  Fig.  5  increases  as  damage  grows.  The  crack  tip  strain  compares  well  with  the 
experimental  observation.  Because  of  a  high  strain  value  at  the  notch  tip,  the  experiment 
could  not  measure  the  notch  tip  strain  directly.  Instead,  the  strain  a  little  distance  away  from 
the  notch  tip  was  measured  and  extrapolated  to  the  notch  tip  for  comparison. 


Fig.  4:  Damage  growth  in  a  uniform  composite  specimen 


Applied  Strain  Crw/m) 

Fig.  5:  Plot  of  induced  notch  tip  strain  versus  applied  strain 


Fig.  6:  Plot  of  damage  at  the  notch  tip  vs.  applied  strain 


Fig.  7:  Normalized  damaged  distance  from  the  notch  tip  as  a  function  of  the  applied  strain 


Figure  6  plots  the  damage  initiation  and  growth  at  the  crack  tip.  The  damage  at  the  notch  tip 
begins  at  the  applied  strain  of  0.012  m/m.  Just  after  damage  initiation,  the  damage  growth 
rate  (defined  as  damage  increase  per  applied  strain)  is  almost  linear  for  a  while  and  the  rate 
decreases  as  the  applied  strain  increases  further. 

Figure  7  shows  the  distance  from  the  notch  tip  where  the  induced  damage  reaches  more  than 
50  percent  of  the  maximum  value.  The  distance  was  normalized  with  respect  to  the  hole 
diameter.  The  plot  indicates  that  the  damage  growth  along  the  minimum  section  of  the 
perforated  specimen  is  not  linear.  That  is,  damage  growth  direction  varies  along  with  the 
applied  strain. {The  major  damage  growth  is  approximately  in  the  45  degree  from  the  applied 
loading  direction  from  the  notch  tip?7 


CONCLUSIONS 

The  micro/macro-approach  was  applied  to  model  and  simulate  damage  initiation  and  growth 
in  a  particulate  composite  which  consisted  of  hard  and  stiff  particles  and  a  soft  and  weak 
rubber-like  matrix.  The  major  damage  mode  called  dewetting  was  modeled  at  the  micro-level 
using  interface  springs  and  a  damage  theory.  The  numerically  predicted  stress-strain  curve  of 
the  composite  and  the  notch  tip  strain  including  damage  compared  well  with  experimental 
data.  The  micro/macro-approach  can  be  used  to  predict  the  strength  of  a  notched  composite 
strength. 
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